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■ Abstract. In this article, we prove convergence of the weakly penalized adaptive discontin- 

uous Galerkin methods. Unlike other works, we derive the contraction property for various 
discontinuous Galerkin methods only assuming the stabilizing parameters are large enough 
to stabilize the method. A central idea in the analysis is to construct an auxiliary solution 
from the discontinuous Galerkin solution by a simple post processing. Based on the auxil- 
iary solution, we define the adaptive algorithm which guides to the convergence of adaptive 
discontinuous Galerkin methods. 



1. Introduction 

The design of adaptive finite element methods based on reliable and efficient a posteriori 
error estimates has been the subject in the past [2j [5j EJ [TTJ [28] . The adaptive finite element 
method consists typically the following successive loops of the sequence 
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Convergence analysis of adaptive finite element methods has been initiated by Dorfler [21] 
who introduced an important marking strategy. Subsequently important theoretical de- 
velopments have been made by many researchers. We refer to [221 EZ] for the work 
on conforming finite element methods, to [Tb] ITS] for the results on mixed finite element 
methods, to [151 E] for the work on nonconforming methods and finally to [2^| 1231 [TP] for 
discontinuous Galerkin methods. On the other hand, the optimality of adaptive finite ele- 
ment method is derived in [9] for two dimensional problems and in [27] for high dimensional 
problems. 

In this article, we focus on the low order adaptive discontinuous Galerkin (DG) methods. 
Karakashian and Pascal [23] were the first to prove contraction properties for the symmetric 
interior penalty Galerkin (SIPG) method. Therein, the authors have proved the contraction 
property for SIPG method under an interior node property. Subsequently the interior node 
property is relaxed independently in the works of [23] and [10]. Moreover the quasi-optimal 
convergence rates are derived in [10]. However, the common issue with the three articles 
[2l| 123] [TO] is that the contraction property is derived assuming the penalty parameters are 
sufficiently large (i.e. larger than what is needed for stability of the method). In this arti- 
cle, we prove contraction properties for various symmetric weakly penalized discontinuous 
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Galerkin methods only assuming that the penalty parameters are large enough to guaran- 
tee stability of the method. For example, in the case of the LDG method the stabilizing 
parameters only have to be positive. This is achieved by a new marking strategy that uses 
an auxiliary solution obtained by post-processing the discontinuous Galerkin solution which 
turns out to be the Crouzeix-Raviart non-conforming approximation (20]. In fact, we bor- 
row the marking strategies (151 [8] that have been developed for non-conforming methods and 
show that these are enough to contract the error of the entire DG approximation. 

The weakly penalized method differ from classic DG methods in the fact that only the 
lower moments of the jumps are penalized on interfaces of the triangulation. For example, 
for piecewise linear elements and the SIPG method the penalty term looks like 

Er /KIM- 

In contrast, in the weakly penalized case, one uses the term 

Ef /n e (R])n.(IM) 

eee h ne Je 

where we use the average of the jump Il e ([[wft]) = ^- / e |iy/i]. Note that this is equivalent to 
using the midpoint rule to evaluate the integrals / so the weakly penalized method 

is cheaper to implement. Weak penalization has been used in weakly over penalized methods 

[I21CE3]. 

We consider the following model problem of finding u G Hq(Q) such that 

(1.1) a(u,v) = (f,v) VveH^Q), 
where 

(1.2) a{w, v) = {Vw, Vt») Vw, v G H^(tt), 

and (•, •) denotes the L 2 (Q) inner product. We assume that Q C M 2 is a bounded domain 
with polygonal boundary dQ and / G L 2 (Q). 

The rest of the article is organized as follows. In Section 121 we introduce the notation 
and preliminary results. In Sectional we recall the DG methods and corresponding stability 
results. In Section HJ we construct an auxiliary solution by averaging the DG solution and 
there in we derive some useful properties and results for the auxiliary solution. Section [5] 
is devoted to the convergence analysis of DG methods. Finally we conclude the article in 
Section [6j 

2. Notation and Preliminaries 
The following notation will be used throughout the article: 

7~h = a face to face, shape regular simplicial triangulations of Q 
T = a triangle of Th hr — diameter of T 
E % h = set of all interior edges of Th 
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E h = set of all boundary edges of Th 

M.t = set of midpoints of the edges of T 
M.\ = set of all midpoints of the edges in £ % h 
A4 b h = set of all midpoints of the edges in £ h 

M h = M l h UM b h 
h e = length of the edge e££/, 
Vh = piecewise (element-wise) gradient 
P m (T) = space of polynomials of degree less than or equal to m > and defined on T. 
The discontinuous finite element space is defined by 

V h = {v h G L 2 (Q) : v h \ T G P 1 (T)}. 
In the analysis below, we need the following Crouzeix-Raviart nonconforming space [20] : 



V CR = {v h eV h : f[[v h jds = VeG^}, 

J e 



and the following vector valued discrete space: 

W h = {w h G L 2 {9) 2 : w h \ T G [P (T)] 2 }. 
Define a broken Sobolev space 

H\n 1 r h ) = {veL 2 (n): Vt = v\ t eh 1 (t) vtgTJ. 

For the DG methods, we require to define jump and mean of discontinuous functions. For 
any e G E\, there are two triangles T + and T_ such that e = dT + fl <9T_. Let be the unit 
normal of e pointing from T_ to T + , and n + = — (cf. Fig. 12.11) . For any v G H l (Q, Th), 
we define the jump and mean of v on e by 

[t>]] = t>_n_ + v + n + , and {{f }} = ~^ v +^' res P ec tively, 

where v± = v\ T . Similarly define for w G H 1 ^, Th) 2 the jump and mean of w on e G £ l h by 

\w\ = u»_ ■ n_ + w + ■ n + , and {{w}} = + w + ), respectively, 

where w± = w\t±- 

For any edge e G there is a triangle T e Th such that e = <9T fl <9f2. Let n e be the unit 
normal of e that points outside T. For any t> G i/ 1 (T), we set on 6 G 

[[f ] = vn e and {{f }} = f , 

and for w e H^T) 2 , 

[wj — w • n e , and {{w}} = w. 
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Figure 2.1. Two neighboring triangles T_ and T + that share the edge e = 
<9T_ n dT + with initial node A and end node B and unit normal n e . The 
orientation of n e = n_ = — n + equals the outer normal of T_, and hence, 
points into T+. 

The discontinuous Galerkin methods use a lifting operator r : L 2 (Sh) 2 — >■ defined by 

(2.1) / r(w) • r = - E / w ■ {{r}} ds Vr G W A , 

and a local analogue r e : L 2 {e) 2 — > Wh defined by 

(2.2) / r e (w) -Tdx = - I w ■ {{r}} ds Vr G W h . 

Jn Je 

Let Il e : L 2 (e) — > R be the L 2 -projection onto constants defined by 
(2.3) 



n e (w) = t- / f ds. 

he Je 

3. Discontinuous Galerkin Methods 

We consider four stable and symmetric weakly penalized discontinuous Galerkin Methods 
namely, the IP (or SIPG method) method [22J [29l [3] , the LDG method pU H], the method 
by Brezzi et al. [H] and the method by Bassi et al. [7] . The original articles have considered 
the formulations using the penalty term with strong jumps, we replace them here with weak 
jumps. The concept of stabilizing the DG formulation by weak jumps was introduced in 



ds 



The bilinear form for the IP method [221 EHJ E] is defined by 
(3.1) A h {w h , v h ) = E J ^ w h ■ V^h dx - JfiVwh 



TeT h 



e££ h 



E / WW ds + E f / n e (Kl)n e (KJ). 



where a > is the stabilizing parameter. 
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The bilinear form for the LDG method [191 H] is defined by 

(3.2) A h (w h , v h )=J2 [ Vw h -Vv h dx-J2 I {{Vw h }}lv h l ds-J^ I {{VMMM ds 

TeT h Jt ee£ h Je ee£ h Je 

+ f r(n e (M))r(n e (M)) + ^f f n.(M)n.(NB). 

The bilinear form for the Brezzi et al. method [H] is defined by 

(3.3) A h {w h , v h )=J2 [ -Vv h dx-J2 [ {{Vw4} M ds-J2[ ds 

TeT h Jt ees h Je eee h Je 

+ / r(n e (W))r(n e (W)) + ^ a / r e (n e (K]]))r e (n e (M)). 

Ja ees h Jn 

The bilinear form for the Bassi et al. method [7] is defined by 

(3.4) A h (w h , v h )=J2 / Vw h -Vv h dx-J2 / ftywhbbhl ds-J2 ii Vv h}}l w hl ds 
' ' Jl ees h ■>' ees h ■'■ 

,(n.(tt«; fc ]l))r.(lI B (K]))- 



Ten 



ee£ h 



a / r R 

n 



The DG method is to find Uh € Vh such that 

(3.5) A h {u h ,v h ) = (f,v h ) Vv h eV h , 

where Ah is any of the bilinear form defined in (I3.ip - (I3.4|) . 

It is proved in JTJ Lemma 1] that the IP method is stable for any a satisfying 

(3.6) a > 4 maxp(SV), 

Ten 

where p(St) is the spectral radius of the local stiffness matrix [5y] mn = (V^A™, VhA n ), where 
A m 's are barycentric coordinates of T. In the Table 13. 1\ we present the condition on a for 
the above DG methods. 



Method 


Condition on a 


IP Method [HI Hi E] 
LDG Method [H ID 
Brezzi et al. [14] 
Bassi et al. [7] 


a satisfies (|3.6p 
a > 
a > 
a > 3 



Table 3.1. Conditions for a 



Define the mesh dependent norm 

(3.7) |M|? ifc = /|v*f^ + ]Tn e (H) 2 wveH\n,r h ). 



Ten ' 



ee£ h 
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The following lemma on the stability of the DG methods (13. 5p is well-known [3]. 
Lemma 3.1. Assume that a satisfies the conditions in Table \3A\ Then, it holds that 

C\\v h \\l h < A h {v h ,v h ) Vv h eV h . 

4. An Auxiliary Solution by Post Processing 

Let M/, G 14 be the solution of any of the DG methods (13.51) . Define an auxiliary solution 
u* h G Vcr by the following: 

{{{u h }}(m e ), if m e e M\ 
0, if m e e M b h . 

In the following lemma, we establish an integral relation for Uh and u* h . 
Lemma 4.1. For any Vh € Vh, it holds that 

(4.2) £ I Vu h ■ Vv h dx-J2 [ {{V^»M ds = [ V < • Vv h ds. 

TeT h ^ T e££ h ^ e TeT h ^ T 

Proof. Using integration by parts, we find 

E / V(%- <) • Vv h dx = / {{Vv h }}[[u h - <]] ds + J^ iVufcHw/, - <}} ds. 

The definition of u* h implies 



E Av^MW-<}}^ = o, 



E fivaKl^ = o. 

This completes the proof. □ 

In the following lemma, we estimate the error between Uh and u* h . 
Lemma 4.2. It holds that 

E (^ll^-^ll^m + llvK-on^^)) <cE n e 



Proof. The proof is an easy consequence of the following estimate: for any Vh € P\(T), it 
holds that 



(4-3) Killer) < C\T\ E Vh(m e ) 2 . 



m e £MT 

□ 



The following identity is useful in our subsequent analysis. 



uAds 
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Lemma 4.3. It holds that 

A h (u h - u* h , u h - u* h ) = (/, u h - u* h ) 

Proof. Using f !3.5j) . we find 

A h (u h - u* h , u h - u* h ) - (/, u h - O = -A h (u* h , u h - u* h ) 

It is remaining to show that Ah{u* hl Uh — u* h ) = 0. Using the fact that II e ( [it£]) = for all 
e G £h, integration by parts and (14. ip . we obtain 

A h (u* h ,u h -u* h ) = E / v < • V K - K) dx - E /{{ V ^}}K 

= E A v <MW-<»^ 

= o. 

This completes the proof. □ 

In the following lemma, we estimate the error u* h — Uh by volume residual. 

Lemma 4.4. Assume that a satisfies the conditions in Table 13 . 11 Then there exists some 
C* > such that 

\\<-n h \\l h <C*J2 h T\\f\\h(Tr 
TeT h 

Proof. Using Lemma 13. 1[ Lemma I4.3[ Cauchy-Schwarz inequality, Lemma 14.21 and Young's 
inequality, we find 

C\\u* h - u h \\ 1;h < A h (u* h - u h ,u* h - u h ) = (f,u* h - u h ) 

/ \ 1/2 / \ 1/2 

< I E h T\\f\\h(T) I E h T 2 W< ~ U h\\l*{T) 

\TeT h J \TeT h 

1/2 / \ 1/2 



< 



< 



c(j2 h T\\nh iT) ) fE n 

7(E^II/II1 2( t))+«-^|IU. 



2 



We complete the proof by choosing e sufficiently small. □ 
Lemma 4.5. T/ie auxiliary solution u* h satisfies 

(V h u* h ,V h v h ) = (f,v h ) \/v h e V CR , 
i. e, u* h is the solution of the classical nonconforming method [20J . 
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Proof. Using ( 13. 5 j) for any Vh G Vc-h, we find 

Then using (14. 2p . we complete the proof. □ 

The following a posteriori error estimator is an easy consequence of Lemma 14.41 and the 
results in [15]: 

Lemma 4.6. Let u anduh be the solutions of (jl.ip and (13. 5p . Letu* h be the auxiliary solution 
defined in (14. ip . JTien /ioWs i/ierf 

|V fc («-tt h )llL»(n)<C7(x;^ll/lli»cr, + E f Ud</ dsf ds 

where d/ds denotes the tangential derivative along the edge e. 
Proof. First using triangle inequality 

||V fc (lt - U h )\\ L 2 {n) < \\V h (u - 0||l2(Q) + \\V h (u* h - Ufc)||i,a(n), 

and then using Lemma 14.41 and the results in fTS] , we complete the proof. □ 



5. Convergence of Adaptive DG Methods 

Let Uh G Vh be the solution of any of the DG method (13.51) and let u* h G Vcr be the 
auxiliary solution constructed in (14. ip . Using Lemma I4.5[ recall that u* h is the solution 
of the Crouzeix-Raviart nonconforming method and Lemma 14.41 implies that there exists a 
positive constant C* > such that 

\K-u h \\l + J2 u ^MY<c*\\hf\\\ 

e££ h 



where hereafter || • \\ h = ||Vj, ■ || and 

" ' U ll T\\J \\L 2 (T)- 



\\hff = hi 11 f " 2 



T£T h 



h 



Let Th be the conforming refinement of Tu obtained by refining the all the marked elements 
in Th that are marked in the step MARK. The functions with h (resp. H) suffix corresponds 
to the mesh Th (resp. Th)- Below, we consider separately two different marking strategies 
that are introduced by Carstensen and Hoppe [15J and by Becker, Mao and Shi [8] and prove 
the error reduction for both the algorithms separately. 

Marking by Carstensen and Hoppe 
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Given the universal constants with < Q, p 2 < 1, the outcome of MARK is a set of edges 
A4 C Sh such that 

(5.1) 6 [ H e ldu* H /dsfds < [ H e ldu* H /dsfds. 

The refined regular triangulation Th from REFINE generated by refining at least all the edges 
in Ai (and possibly further edges to avoid hanging nodes) with the new mesh-size h < H is 
supposed to satisfy 

P2\\Hf\\l Hn) < \\hf\\l Hny 

The results in Carstensen and Hoppe [15] imply that there exists < p± < 1 and C\ > 
such that 

(5-2) ||u - u* h f h < Pl \\u - u* H \\ 2 H + Ci\\Hf\\ 2 , 

(5-3) ||^|| 2 <p 2 ||tf/|| 2 . 

In the following theorem, we derive the contraction property for the adaptive DG methods 
(13. 5p using the marking strategy by Carstensen and Hoppe (15.11) . 

Theorem 5.1. Let the marking be done by (15.11) . Then there exists 7 > and < p* < 1 
such that 

h - uhU + l\\hf\\ 2 < P* {h - u H \\l + l\\Hft) • 
Proof. Let e > 0. Using triangle inequality and Young's inequality, we find 

\\u - u h \\l + 7 ||/>/|| 2 < (1 + e)\\u - ul\\l + (1 + l/e)|K - u* h \\ 2 h + ^hff. 



Using Lemma [4.41 we obtain 

\\u - u h \\l + 7 ||/*/|| 2 < (1 + e)||u - + (C*(l + l/e) + 7) \\hf\\ 2 . 
Using the error reduction for u* h (15. 2p . we find 

\\u - u h \\l + 7IMI 2 < (1 + e) Pl (\\u -u* H \\ 2 H + C.WHft) + (C*(l + l/e) + 7) \\hf\\ 2 . 
Again using triangle inequality, Young's inequality and lemma I4.4[ we find 

\\u - u h \\l + j\\hf\\ 2 < (1 + e) Pl ((l + e)\\u - u H \\ 2 H + C*(l + l/e)\\Hf\\ 2 + C.WH f\\ 2 ) 

+ (C*(l + l/e) + 7 )||/ i /|| 2 . 

Therefore 

\\u - u h \\ 2 h + 7 |M| 2 < (1 + e) 2 Pl \\u - u H \\ 2 H + (1 + e)pi (C*(l + l/e) + C x ) \\Hf\\ 2 
+ (C*(l + l/e) +1 )\\hf\\ 2 . 
We now use (15.31) and find 

\\u - u h \\ 2 h + 7 ||/>/|| 2 < (1 + e) 2 P i\\u - u H \\ 2 H + (1 + e) Pl (C*(l + l/e) + C x ) \\Hf\\ 2 

+ (C*(l + l/e) +1 )p 2 \\Hf\\ 2 . 
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By simplifying 

||u - u h \\l + i\\hf\\ 2 < (1 + ef Pl \\u - u h \\h 

+ [(1 + e)pi (C*(l + 1/e) + d) + (C*(l + l/e))p 2 + 7 p 2 ] 1 1#/|| 2 . 

Choosing 7 sufficiently large such that 

[(1 + e) Pl (C*{1 + 1/e) + d) + (C*(l + l/e))p 2 + 7P2] < 7(1 + p 2 )/2. 

That is by choosing 7 such that 

2 [(1 + e)pi (<T(1 + 1/e) + d) + (<T(1 + l/e))p 2 ] /(l - p 2 ) < 7, 

we complete the proof by p* = min{(l + e) 2 pi, (l + p 2 )/2} and e such that (l + e) 2 pi < 1. □ 

Marking by Becker, Mao and Shi [Hj: Choose the parameters < 6, a < 1 and 7. Then the 
out come of MARK step is the set of edges or elements according to the following: 
If \\Hf\\ 2 L 2( n j < lYlie^ f e H e ^du* H /ds]\ 2 ds, then mark a subset M. C E\ with minimal 
cardinality such that 

(5.4) 6 [ H e ldu* H /dsfds < [ H e ldu* H /dsfds, 

else find a subset Ai\ C Th with minimal cardinality such that 

^tII/IIl^t)- 



(5.5) a H 2 T \\f\\l^T)< £ 



TgTh T£Mi 

The result [8] by Becker, Mao and Shi is that there exist < p < 1 and /3* > such that 

(5.6) || U - <|| 2 + p\\hf\\ 2 < P(h - u* H \\ 2 H + /3||tf/|| 2 ), 

for all such that 0>0*. 

Below, we prove the convergence of adaptive DG methods under the Becker, Mao and Shi 
marking fl5.4l) - fl5.5p . 

Theorem 5.2. Suppose that the marking is done by (15.4p - (15.5H . Then there exists 7 > 
and < p* < 1 such that 

\\n - «fc||J + 7ll^/l| 2 < P* (\\u - u H \\ 2 H + 7 ||# Z|| 2 ) • 

Proof. Let e > 0. Using triangle inequality and Young's inequality, we find 

\\u - u h \\l + l\\hf\\ 2 < (1 + e)\\u - <|| 2 + (1 + l/e)\\u h - <|| 2 + 7 ||/i/|| 2 . 



Using Lemma I4.4[ we obtain 

\W - u h \\ 2 h + 71^/H 2 < (1 + e)||u - ul\\l + (C*(l + 1/e) + 7) \\hf\\ 2 , 
or equivalently 

\\u ~ u h \\l + 7||^/|| 2 < (1 + e) (\\u - <|| 2 + (C*(l + 1/e) + 7 )/(l + e)||/>/|| 2 ) . 
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Assume that 7 is sufficiently large such that 

(5.7) (C*(l + l/e)+ T )/(l + e)=:/9>/3*. 

Then, using ( 15. 6 p we arrive at 

\\u- Uh \\i+ ^hff < (1 + e) P (\\u -u* H f H + mn 2 )- 



Again using triangle inequality and Young's inequality and Lemma I4.4[ we find 

\\u ~ u h \\l + 7 ||/*/|| 2 < (1 + e)p [{I + e)\\u - u H \\% + C*(l + l/e)\\Hff + /3||tf/|| 2 ] . 
Therefore 

||« - u h \\\ + 7 ||/i/|| 2 < (1 + efp\\u - u H \\ 2 H + (1 + e)p(C*(l + l/e) + 0) \\Hf\\ 2 . 

First note that we can choose e sufficiently small such that (1 + e) 2 p < 1. Then the proof 
will be completed if we can show that there is some < P2 < 1 such that 

(5.8) (l + e)p(C*(l + l/e) + /3) < p 27 , 

with p* = min{(l + e) 2 p, p 2 }. 
Using fl52D in (l5T8|) . 

% , , , (l + e)pC*(l + l/e) 
1 + e ) P C* 1 + l/e + i '-j— y - '-L + p 7 < p 27 , 

(1 + 6) 



equivalent ly 



:i+6)P c-( 1+1 / e)+ ' i+ y;(; + ^ <(p 2 -ph, 



The proof is completed by choosing p 2 such that p < p 2 < 1 and 7 sufficiently large. □ 
Remark 5.3. Using triangle inequality and Lemma |4.4( we find 

||« - Uh\\h < \\u - U* h \\ + \\u* h - Uh\\h 

<\\u-ul\\+C\\hf\\. 

Therefore the adaptive DG methods converge at least at the rate of adaptive nonconforming 
method. 

6. Conclusions 

In this article, we have proved the contraction property for various symmetric discon- 
tinuous Galerkin (DG) methods. Unlike in the existing works for strongly penalized DG 
methods, we prove the convergence of weakly penalized adaptive DG methods without fur- 
ther assuming the stabilizing parameter is larger than what is required for stability. Although 
the analysis in this article is restricted to the lowest order case, we hope that similar ideas 
may be used in higher order cases. We remark that the convergence analysis of adaptive 
DG methods using strong jumps is still open when the stabilizing parameter is chosen just 
according to the stability of the method. 
Acknowledgments 
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